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§ ABSTRACT 
(N 

I We have determined the ratios of total to selective extinction in the near- 

infrared bands (J, H, Ks) toward the Galactic center from the observations of 
the region | / |< 2?0 and 0?5 <| 6 |< 1?0 with the IRSF telescope and the SIRIUS 
^ _ camera. Using the positions of red clump stars in color-magnitude diagrams as 

' a tracer of the extinction and reddening, we determine the average of the ratios 



of total to selective extinction to be Ak^/ Eh-Ks = 1-44 ± 0.01, Ak^/ Ej^Ks = 
0.494 ± 0.006, and Ah/Ej_h = 1-42 ± 0.02, which are significantly smaller than 
^ I those obtained in previous studies. From these ratios, we estimate that Aj : 

O ; Ah : Aks = 1 ■■ 0.573 ± 0.009 : 0.331 ± 0.004 and Ej_h/Eh-Ks = 1-72 ± 0.04, 

and we find that the power law Ax oc ;\-i-99±o o2 jg good approximation over 
these wavelengths. Moreover, we find a small variation in Ak^/ Eh-Ks across 



O 

in ' our survey. This suggests that the infrared extinction law changes from one line 



O ■ 

^ ' of sight to another, and the so-called "universality" does not necessarily hold in 

the infrared wavelengths. 
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1. Introduction 

The wavelength dependence of interstellar extinction is not only important information 
for understanding the nature of interstellar dust grains but an essential ingredient in recov- 
ering the intrinsic properties of reddened objects. Change in the interstellar extinction law 
affects many published conclusions based on the observations of heavily reddened objects, 
such as photometric distances, the selection of young stellar objects in color-color diagrams, 
and the age determination of clusters from comparison with isochrones. Since the Galactic 
center (GC) region is not observable at visible wavelengths due to large interstellar extinc- 
tion, observations and a precise extinction law at infrared wavelengths are required for the 
study of objects in the GG. However, accurate determination of the absolute extinction ^4;^ 
is not easy in the infrared. 

The most widely used technique to obtain the interstellar extinction law is the color- 
difference (GD) method, in which the extinction curves are determined in the form of the 
ratios of color excesses -Ea-Ai/-^A2-Ai as a function of A^^. Here Eb-v (^i = ^ and A2 = 
B) is customarily used as a standard normalization. The absolute extinction A\ can be 
derived from the ratios of color excesses, only if the ratio of total to selective extinction 
Rv — Ay /Eb-v is evaluated. While the ratio of color excesses is easily observed and thus 
determined precisely, Ry cannot be evaluated directly from observed data; Ry can only be 
deduced by the extrapolation of the extinction curve to = with reference to some model 
(e.g., van de Hulst 1946) for its behavior at wavelengths beyond the range. The intercept of 
the extinction curve in the E\-v / Eb-v versus A~^ diagram is —Ry because 

hm 1^ = ^ = -Ry. (1) 

A 1— >0 Eb-V Eb-V 

In this diagram, then, the total extinction (normalized by Eb-v) at each wavelength A is 
obtained as the vertical distance Ax/ Eb-v of each data point from the level of —Ry through 
the relation 

^^-^ {-Rv). (2) 



Eb-v Eb-v 

Unfortunately, the extrapolation has uncertainty due to (1) the contamination by dust emis- 
sion and (2) the possible existence of neutral extinction by grains much larger than the 
wavelength of observation. The total extinction value A\ is derived as the difference be- 
tween the two terms on the right-hand side^ of eq. [2] . Since the value of eq. [2] is small in 



^Note that Ex-v is negative for the wavelengths A longer than V. 
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the infrared, even a small change in Ry can have significant infiuence on the value. Therefore, 
the accurate determination of is difficult in the infrared. 

Red clump (RC) stars in the Galactic bulge have recently been used to study the ex- 
tinction law (Wozniak & Stanek 1996), resulting in a precise determination of the ratio of 
total to selective extinction Rx (Udalski 2003; Sumi 2004). RC stars are the equivalent of 
the horizontal-branch stars for a mctal-rich population, and they have narrow distributions 
in luminosity and color with weak dependence on metallicity. In addition, distances to RC 
stars in the Calactic bulge can be regarded as the same. They thus occupy a distinct region 
in the color- magnitude diagram (CMD), and can be used as standard candles at an equal 
distance; they are located in a straight line with a slope of R\ = Ax/E^>_-^, in accordance 
with the variable extinction to each of them. We refer to this as the "RC method" . The 
RC method is based on the "variable extinction method" (Krelowski & Papaj 1993) or the 
"cluster method" (Mihalas & Binney 1981). 

In the cluster method, we expect the difference between the observed and absolute 
magnitudes to be the same for every cluster star except for differences produced by variations 
in the amount of interstellar extinction along the line of sight to each star (Mihalas & Binney 
1981). Therefore, photometry and the determination of the spectral type of each star is 

necessary to derive the slope R\ = Ax/E^_y Here the amount of interstellar extinction 
is different, but the interstellar extinction law, and furthermore the property of interstellar 
dust, is assumed to be identical between the lines of sight. In the real world, we should note 
that the extinction Ax might be nothing but an average of possibly variant extinction due 
to variant dust grains integrated over the whole line of sight. 

In the RC method, while this assumption is similarly applied, we select RC star candi- 
dates from the CMD and do not have to perform spectroscopy. We can obtain Rx directly 
from the observed data, and therefore Rx can be much more reliable than that obtained by 
the extrapolation of the extinction curve in the CD method. The reliable Rx provides the 
absolute extinction Ax which is exactly what we want to know. While visible observations of 
RC stars are restricted to a few windows of low extinction, we can observe highly reddened 
RC stars at fields very close to the GC in the infrared. Thus, we can use the RC stars near 
the GC to determine the infrared extinction law very accurately and compare the results 
with the past determinations for the GC and other lines of sight. 

In this paper, we have determined the ratios of total to selective extinction Axg/ Eh-Ks-i 
AksI Ej_Ks-i and Ah/Ej_h, and the ratio of absolute extinction Aj : Ah : A^g, by adopting 
the RC method to investigate the extinction law toward the Galactic center. 
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2. Observations 

Our observations were obtained in 2002 March— July, and 2003 April— August using 
the near-infrared camera SIRIUS (Simultaneous Infrared Imager for Unbiased Survey; Na- 
gashima et al. 1999; Nagayama et al. 2003) on the IRSF (Infrared Survey FaciUty) telescope. 
IRSF is a 1.4m telescope constructed and operated by Nagoya University and SAAO (South 
African Astronomical Observatory) at Sutherland, South Africa. The SIRIUS camera can 
provide J(1.25//m), i?(1.63/im), and Ks{2.1Aiim) images simultaneously, with a field of view 
of 7'.7 X 7'.7 and a pixel scale of ^.'45. 

Over the period 2002-2003, about 800 x 3 (J, H, Kg) images were obtained over | / |< 2?0 
and I 6 |< 1?0 (see Figl). Our observations were obtained only on photometric nights, and 
the typical seeing was 1".2 FWHM in the J band. A single image comprises 10 dithered 5 s 
exposures. 




6=0° 



(Ks vs. H-Ks) 



Fig. 1. — Observed area and fields used for data analysis. Each square consists of nine 
SIRIUS fields. Only the hatched regions were used for the data analysis in this paper. The 
white squares are the regions where the magnitude and color of RC stars were not reliably 
determined due to large extinction. 
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3. Reduction and Analysis 

Data reduction was carried out with the IRAF (Imaging Reduction and Analysis Fa- 
ciUty)^ software package. Images were pre-reduced following the standard procedures of 
near-infrared arrays (dark frame subtraction, flat-fielding, and sky subtraction). Photom- 
etry, including point-spread function (PSF) fitting, was carried out with the DAOPHOT 
package (Stetson 1987). We used the daof ind task to identify point sources, and the sources 
were then input for PSF-fitting photometry to the allstar task. About 20 sources were 
used to construct the PSF for each image. 

Each image was calibrated with the standard star #9172 (Persson et al. 1998), which 
was observed every hour in 2002 and every half-hour in 2003. We assumed that #9172 has 
J = 12.48, H = 12.12, and Kg = 12.03 in the IRSF/SIRIUS system. The average of the 
zero-point uncertainties was about 0.03 mag in the three bands. The averages of the lOcr 
limiting magnitudes were J — 17.1, H — 16.6, and Kg — 15.6. 

To measure the ratios of total to selective extinction Aks/Eh-Ks, -^KsI Ej-Ks^ and 
Ah/Ej^h, we selected the bulge RC stars, which constitute a compact and well-defined 
clump in a CMD. Their mean magnitude and color are thus good references, and we apply 
the method developed in the V and / bands by Udalski (2003) and Sumi (2004). 

As a first step we divide each field into four sub-fields of ~ 4' x 4' on the sky. Then we 
construct Ks versus H — Ks, Ks versus J — Ks, and H versus J — H CMDs. Second, we 
extract stars in the region of CMDs dominated by RC stars (the rectangular region of the 
CMD in Fig. 2), and the extracted stars are used to make magnitude (luminosity function: 
Fig. 2, upper right) and color (Fig. 2, lower left) histograms. These histograms have clear 
peaks that can be fitted with a Gaussian function (thick curves in Fig. 2 histograms). 

Due to highly nonuniform interstellar extinction over the area surveyed, the RC peaks 
in CMDs shift from one sight line to another. The peaks shift in the range 13.4 < Ks ^ 14.6 
and 0.4 < H — Ks ^ 1.2, and therefore we have to shift the region to extract RC stars 
from sub-field to sub-field. Since the mean J, H, and Ks magnitudes of RC stars become 
too faint in highly reddened fields, estimates of the peak magnitudes and the colors of RC 
stars can be unreliable in such fields. To avoid this problem, we use only the sub-fields in 
which the peak magnitude of RC stars is more than 1 mag brighter than the lOcr limiting 
magnitudes (Fig.l, hatched squares). In addition, we confirmed the completeness to be 85% 



^IRAF is distributed by tlic National Optical Astronomy Observatory, which is operated by the Associa- 
tion of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 
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at Ks = 15 by adding artificial stars into the most crowded image. 




H-Ks [mag] 



Fig. 2. — Upper left: Sample CMD centered at 17h 46m 10.0s, -27d 13m 48.1s. Stars in the 
rectangle on the CMD are used to estimate the RC peak of color and magnitude. Upper 
right: Luminosity function of selected region in CMD. Lower: H — Ks color histogram of 
selected region in CMD. Mean color and magnitude are obtained by fitting with Gaussian 
function. 
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4. Results 

Fig. 3 shows the location of the RC magnitude and color peaks in the Ks versus H — Ks 
{upper left), Ks versus J — Ks {upper right), and H versus J — H {lower) CMDs. Error bars 
include uncertainties in the RC peak and the photometric calibration. Solid lines are least- 
squares fits to the data. If the variations of the RC magnitudes and colors are due to the 
interstellar extinction, the ratios Aks/ Eh-Ks-i -^Ks/Ej-Ks, and Ah/Ej^h can be directly 
determined by the slopes of the fitting hnes. Thus, we obtain Aks/Eh-Ks — 1-44 ± 0.01, 
Aks/Ej-Ks = 0.494 ± 0.006, and Ah/Ej^h = 1-42 ± 0.02. 

AksI Eh-Ksi ^KsI Ej-Ks-i and Ah/Ej^h provide us with the ratios of absolute extinc- 
tion, Aks/Ah, Aks/Aj, and Ah/Aj, by which we obtain Aj : Ah : Aks = 1 : 0.573±0.009 : 
0.331 ± 0.004. Here Ah/Aj = 0.573 is the weighted mean of that derived by Ah/Ej_h and 
that derived by Aks/Eh-Ks ^Ks/Ej-Ks [i-G-, from {Aks/Aj)/{Aks/Ah)]- Using the 
ratios Ah/Ej^h = 1-42 ± 0.02 and Ah / Eh-Ks = 2.43 ± 0.01, in turn, we can determine 
the most frequently quoted "near- infrared extinction law", the ratio Ej_h / Eh-Ks to be 
1.71 ± 0.02. When we plot the observed peaks of RC stars in the J — H versus H — Ks 
color-color diagram, the slope of peak distribution Ej^h/ Eh-Ks = 1-72 ± 0.04 (Fig 4) is 
consistent with this ratio. 

To examine the variation of the extinction law we divide the survey region into four 
sub-regions, N+ (+2° > / > 0°,+l° > 6 > 0°), S+ (+2° > / > 0°,0° > 6 > -1°), N- 
{0° > I > -2°,+l° > b > 0°), and S- (0° > / > -2°,0° > b > -1°). The CMDs of the 
divided regions are shown in Fig. 5. The Ak^/Eh-Ks ratios for N+, S+, N— , and S— are 
1.46 ± 0.02, 1.45 ± 0.02, 1.34 ± 0.02, and 1.51 ± 0.03, respectively. We can see a small but 
clear difference from one line of sight to another. As discussed in §5.3, the extinction curve 
can be approximated by a power-law A\ oc A~". Then the power-law exponents a for N-|-, 
S+, N— , and S— are 1.96, 1.97, 2.09, and 1.91, respectively. The low numbers of data points 
in CMDs of Ks versus J — Ks and H versus J — H create uncertainties that are too large 
to detect variations of similar size to those found in Axg/ Eh-Ks- 
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Fig. 3. — Location of RC peaks in Kg vs. H — Kg {upper left), Ks vs. J — Ks {upper right), 
and if vs. J — H {lower) CMDs. The solid lines are the least-squares fits to the data. 
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0.3 0.4 0.5 0.6 0.7 0.8 

H-Ks [mag] 



Fig. 4. — Location of the RC peaks in the color-color diagram { J — H vs. H — Ks )■ The 
solid hne is the least-squares fit to the data. 



-10- 




Fig. 5. — Variation of Ak^/ Eh-Ks- The solid lines are the least squares fits to the data. 
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5. Discussion 

5.1. Possible Errors in R 

An advantage of the RC method is the large number of RC stars, which makes their 
mean magnitude determination straightforward and precise. RC stars are also numerous 
in the solar neighborhood. Alves (2000) derived the average K band luminosity of 238 
Hipparcos RC stars, Mk — —1.61 ± 0.01, and their precise calibration is possible. 

Population effects on the intrinsic magnitude and color of RC stars are the main error 
sources of the method. It is known that their brightness depends on their metallicity; how- 
ever, the dependence is very weak (< 0.05 mag / [M/H]; Salaris & Girardi 2002). To change 
the slope that we obtain, the metallicity distribution of RC stars in the bulge would have to 
be highly correlated with the extinction. Moreover, no metallicity gradient has been found 
in I / |< 4° (Frogel et al. 1999; Ramirez et al. 2000). 

Another error source is the variation of distance to RC stars, changing their mean magni- 
tude by ~ 0.05 mag along the Galactic longitude in the current observation range (Nishiyama 
et al. 2005). This difference can change the slope Aks/Eh-Ks of 1-44 [/^Ks/^{H - Ks) = 
1.37/0.95, see the Kg vs. H - Kg CMD in Fig. 3] to 1.49 (=1.42/0.95). The change is only 
0.05, even if the distance is completely correlated with the extinction, which is very unlikely. 
Therefore we conclude that the variation of the metaUicity and distance does not affect our 
results. 

The change in the extinction law, as shown in Fig. 5, suggests that within a single field, 
there may also be similar variations that would produce a selection bias in the results. To 
estimate the size of such bias, we divide the data points plotted in the Ks versus H — Ks 
CMD into two groups: those having a counterpart in the Ks versus J — Ks CMD (bluer 
sources) and those without one (redder sources). The former has a color distribution 0.4 < 
H - Ks < 0.7 and a slope Ak JEh-Ks = 1-46, and the latter has 0.5 < H - Ks < 1.2 and 
Ak<^/ Eh-Ks = 1.43. This means that the selection bias due to the extinction, which probably 
reflects the distribution of sources along the lines of sight, produces an uncertainty of only 
a few per cent in our results. This small uncertainty is also confirmed by the consistency 
of Ej_h/Eh-Ks ~ 1-7 derived by the distribution of RC stars in the color-color {Ej_h vs. 
Eh-Ks) diagram (Fig. 4) and derived by the product of Ah/Ej_h and Ah/Eh-Ks- Note 
that the sources whose J magnitudes have been determined form a small sub-group. 

The IRSF/SIRIUS system is similar to the MKO system (Tokunaga et al. 2002), and the 
effective wavelengths following equation (A3) of Tokunaga & Vacca (2005) are calculated for 
a typical bulge RC star when it suffers extinction of Ak^ ~ 0.4—1.6 (the range of extinction in 
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Fig. 3). We have employed a model of Kurucz (1993) with the parameters of Teflf = 4750K, 
[M/H] = —0.1 dex, and logg = 2.0, and a power-law extinction curve (see §5.3). The 
effective wavelengths change by ~ 0.01/xm to the longer wavelengths as the star suffers more 
extinction, and the mean effective wavelengths of the J, H, and Kg filters are 1.25, 1.64, and 
2.14 fim, respectively. The changes are small, and these wavelengths are not very sensitive to 
the source spectral energy distribution, and we thus conclude that the effective wavelengths 
can be uniquely defined for our study. 



5.2. Compcirison with Previous CD Method Studies 

Previous results in the literature are listed for comparison in Table 1. Our results show 
a clear difference from those obtained in previous work. The ratios of total to selective 
extinction determined in this study are smaller, which corresponds to faster decrease in the 
absolute extinction when the wavelength increases. Therefore, the resulting power-law 
index a is larger in this study. 

The theoretical curve No. 15 of van de Hulst (1946, hereafter vdH) based on the Mie 
scattering theory has been used to estimate the wavelength dependence and is still favored by 
several authors (e.g.. Glass 1999; Jiang et al. 2003). Our results are the closest in agreement 
to vdH. 

Rieke & Lebofsky (1985, hereafter R&L) determined the near-infrared extinction law 
by applying the CD method to their observations of five supergiants near the GC and o Sco. 
Their results have quite often been referred to as the standard extinction law, especially 
toward the GC, so let us examine them. Their observations were made in the K band, 
not Ks, but we can use results of Cardelli et al. (1989, hereafter CCM89) derived by the 
analytic formula of the average extinction law for 0.9/im < A < 3.3/im by fitting the data 
of R&L with a power law. Since the CD method requires the ratio of total to selective 
extinction, R&L set Ey_K/ Eb~v = 2.744 and assumed Ry = Av/Eb-v = 3.09 for the GC 
sources to evaluate Ax/Eb^v s'-iid Ax/ Ay- While their assumption of Ry is reasonable, a 
small change in Ry results in a large difference in Ax/Eb-v at near-infrared wavelengths 
(see eq. [2]). If Ry toward the GC is Ry — 3.00, which is smaller by only 0.09 than that 
in R&L, we obtain the ratio Aj/Ay : Ah/ Ay : Ax/Ay = 0.260 : 0.150 : 0.085 with their 
Ex-y / Eb-v- These ratios are in fact similar to the ratios shown by vdH, Aj/Ay : Ah /Ay : 
Ak/Av = 0.245 : 0.142 : 0.081. Although the lower limit of Ry in R&L was determined 
by the extinction at L. M, 8/im, and 13/im, the observations at these wavelengths have 
relatively large uncertainties. We also note that the extinction law that they derived in the 
wavelength range L, M, and 8/xm is inconsistent with those obtained toward the GC, in 
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Tablc 1: The wavelength dependence of the interstellar extinction. 
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References. — (1) van de Hulst (1946); (2) Rieke & Lebofsky (1985); (3) CardeUi et al. (1989); (4) He et 
al. (1995); (5) Indebetouw et al. (2005) 

particular by spectroscopy with the Infrared Space Observatory (ISO) for A > 2.4:nm (Lutz 
et al. 1996; Lutz 1999, see §5.3). Therefore, although the inconsistency might be the result 
of spatial variation of the extinction law examined in the previous section, it is also quite 
possible that the inconsistency was caused by the uncertainty of Ry in the CD method. 

Indebetouw et al. (2005) derived Ax/Axg by the CD method, using data of the Spitzer 
Space Telescope (SST) and Two Micron All Sky Survey (2MASS) along two lines of sight 
{I — 42° and 284°) in the Galactic plane. They observed RC stars and determined Ax/Ak^ 
for 1.25//m < A < 8.0/im. They measured color excess ratios Ex-KsI Ej-Ks by fitting the 
locus of RC stars in color-color diagrams. To derive Ax/ A^g from Ex-KsI Ej-Ks^ they first 
determined Aj/Axg by fitting the RC locus in a CMD. Since the fields of / = 42° and 284° are 
dominated by the disk stars, which are located at different distances from us, they assumed 
that the extinction per unit distance is constant to fit the RC locus. Their results are also 
different from ours. We notice, however, that the intercept of Ex-KsI Ej-k at A~^ = in 
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their Fig. 5 seems to be larger (less negative) than their result obtained by fitting the RC 
locus in CMD —Axg/ Ej_k = —0.67, but rather close to our result —0.494. Although the 
data of the SST have the potential to establish the extinction law for 3.55/im < A < 8.0/im, 
the uncertainty of Aj/A^s affects Ax/A^g at these wavelengths, and thus the determination 
of Aj/Aks should be approached cautiously. 
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Fig. 6. — Comparison of extinction ratio (Ax/Axg) toward the Galactic center (this work 
[filled circles] and Lutz et al. 1996 [open circles]), and other regions. The solid line is found 
from the relative intensities of the H2 lines toward the emission peak of Orion OMC-1 (eq. 
2 of Rosenthal et al. 2000). The average of three fields in the Galactic plane (Indebetouw 
et al. 2005) is also shown {open squares). 
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5.3. Power-Law Extinction in the Infrared 

It is frequently remarked that the continuous extinction curve between ~ 0.9 and ~ 5fim 
can be approximated by a power-law A\ oc A~" (Draine 2003). However, the power-law fitting 
has been based on the CD method, and is not free from the offset errors of eq. [1] . In contrast, 
without any a priori assumptions about the extrapolation to the longer wavelengths, our new 
results are fitted very well by the function with a — 1.99 ± 0.02 (see the three data points 
that are very well aligned in Fig. 6). Here we have used the mean effective wavelengths of 
the J, H, and Ks filters determined in §5.1. This good fit presents, for the first time, direct 
evidence for the power-law approximation in the wavelength range of 1.2— 2.2/xm. 

This corroborates the method of determining the absolute extinction on the assumption 
of a power law in the near-infrared without resorting to the CD method. This wavelength 
range has no absorption features so far identified in either the diffuse or molecular interstellar 
medium. Therefore, we proceed with the power-law approximation, which is expressed by a 
straight line in a logarithmic graph like Fig. 6, and compare the index a in the measurements 
in the literature. This also mitigates the difference between various photometric systems used 
in the near-infrared, whose problems have been pointed out by Kenyon et al. (1998). 

He et al. (1995) made photometry of 154 obscured OB stars in the southern Milky Way. 
They basically employed the CD method but derived the absolute extinction also, on the 
assumption that the interstellar extinction follows a power-law. Although their observation 
shows considerable scatter (as suggested by their Fig. 9; also see their Fig. 10), the average 
extinction is clearly different from our results (Table 1). 

Also on the assumption of power-law extinction, Moore et al. (2005) determined the 
power-law index a toward a set of nine ultracompact H II regions and two planetary nebulae 
with 0.27 < r(Br7) < 4.7. Their observations of hydrogen recombination fines at 1.0//m 
< A < 2.2//m showed that the power-law index a varies from 1.1 to 2.0 and a clear tendency 
that a is smaller for higher extinction regions. They suggest that the fiatter near-infrared 
extinction laws might be the result of grain growth in the regions of higher extinction. 
The steep extinction law determined toward the GC in this study, then, might mean that 
the interstellar medium in its line of sight is largely diffuse. A steep index of 2 in the 
polarization, which is supposed to indicate the difference of two orthogonal extinctions, of 
GC objects is also reported by Nagata et al. (1994). Although it has been suggested that 
interstellar extinction remains relatively invariant ( e.g., CCM89 ) in the near-infrared region 
( ^ ^ 0.7/im ), these studies indicate clear variation in the extinction law. 

At longer wavelengths, two sets of extinction determination from the emission line ratios 
using the Short Wavelength Spectrometer on board ISO have been reported toward the GC 
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(Lutz et al. 1996; Lutz 1999) and Orion (Rosenthal et al. 2000). Although these two studies 
have assumed the absolute extinction values at some wavelengths and therefore are vulnerable 
to offset errors, the extinction law toward the GC seems to be clearly different from that of 
Rosenthal et al. (2000), especially in the range 2.5/im ^ A < 8fim (Fig. 6). In the GC Lutz 
(1999) points out that the extinction becomes flatter at A > 2.5/im and lacks the pronounced 
minimum near 7/im expected for the standard grain models; this flatter slope or addition 
of another extinction component is also observed in the polarimetry of the GC objects at 
2.8 — 4.2//m (Nagata et al. 1994). Our extinction law, which has steeper exponents in the 
wavelength range of 1.2 — 2.2//m, and other studies discussed above lead to the conclusion 
that the extinction law varies from one line of sight to another even in the near-infrared 
wavelength range. 



6. Conclusion 

The ratios of total to selective extinction Aks/Eh-Ks — 1-44 ± 0.01, Ak^/Ej^Ks — 
0.494 ± 0.006, and Ah/Ej^h = 1.42 ± 0.02 have been directly obtained toward the Galactic 
center from the observation of bulge red clump stars. Then the ratio of absolute extinction 
Aj : Ah '■ A^^ = 1 : 0.573 ±0.009 : 0.331 ±0.004 is obtained, and the power law approxima- 
tion Ax oc A^^'^^ is shown to be good enough in the wavelength range of 1.2— 2.2/im. These 
values are clearly different from those obtained in previous studies toward the Galactic cen- 
ter and other lines of sight. The previous extrapolation procedures might have caused the 
difference, although the observed color excesses in the literature can be consistent with the 
current results. Furthermore, small variations in the ratio of total to selective extinction do 
exist among the sub- regions in the current study. This indicates that, contrary to the widely 
accepted ideas, the extinction law is not "universal" even in the infrared, and therefore one 
should be extremely cautious about applying reddening corrections. 
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